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1 Introduction 



As has been shown by Mueller 0, the high energy scattering in QCD can be treated in the most 
economical way in terms of color dipole degrees of freedom. In this approach, one considers a fast 
moving particle, as a system of colorless dipoles. The wavefunction of this system of dipoles can be 
found from the QCD generating functional pp. As was shown in Ref. this functional obeys a 
linear functional evolution equation (see also Ref. ji]). This linear functional evolution equation 
was derived in large Nc approximation. For a small projectile, for which we can neglect nonlinear 
effects associated with high dipole densities in its wavefunction, the functional evolution was shown 
to reproduce the dipole version of the Balitsky hierarchy [S] for the scattering amplitude. The 
latter reduces to the Balitsky-Kovchegov (BK) equation El if in addition we assume that the 
low energy dipole interaction with the target has no correlations. 

Though the BK equation has been widely used for phenomenology 0, it is clear starting from 
the first papers on non-linear collective effects at high energies (see Refs. jHUHlE]) that simple 
non-linear evolution equations of the BK type could be correct only in a very limited kinematic 
range. At present there are several approaches allowing one to go beyond the BK equation. 
Balitsky ^ has developed a Wilson line approach which allows the incorporation of both target 
correlations and Nc corrections. This method is equivalent to the effective Lagrangian approach 
describing the Color Glass Condensate and its derivative JIMWLK equation [TT] . 

The methods above describe high energy evolution in a highly asymmetric manner: either 
the projectile or target is always considered as a small perturbative probe. Thus, the constructed 
evolution contains a one way parton shower, in apparent violation of the t channel unitarity ^2] • Is 
is thus challenging to attempt to restore the symmetry of the QCD evolution. If we are successful, 
we would be able to confront problems of high energy scattering of hadrons or heavy nuclei in a 
reliable manner. 

Several steps have been made recently in attempt to formulate a symmetric evolution. Braun 
[T^ used the QCD triple Pomeron vertex ^] both for Pomeron splittings and mergings. lancu 
and Mueller suggested a high energy factorization [T21 CHI UHl, while a statistical approach to 
high energy scattering was proposed in Ref. jT7j. In its turn, Balitsky in Ref. considered a 
symmetric scattering of two shock waves. Another technique is due to Lipatov (e.g. [19j) who 
built an effective theory for reggeized gluons. Unfortunately, a relation between Lipatov' s theory 
and the other methods is so far not clear. We want also to restore a symmetry without loosing the 
probabilistic interpretation of our results which leads to a simplest physical picture of the process 
and a direct application for experimental observations. 

A symmetric way to describe high energy QCD does exist: it is the Reggeon technique [SI El 
El 1201 based on interacting BFKL Pomerons j^Ij. Many elements of this technique are known 
(see Refs. and the main problem is to sum all Reggeon graphs. Past experience in 

summing Reggeon diagrams does not look encouraging. However, a remarkable breakthrough 
was achieved in the last days of the Pomeron approach to strong interactions: it turns out that 
the Reggeon calculus can be re-written in a probabilistic language. It was possible to formulate 
equations for probabilities to find a definite number of Pomerons at fixed rapidity |2S1 
The goal of this paper is to re-write the reggeon calculus of the BFKL Pomerons, by extending 
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our linear functional approach based on dipole generating functional [21 E|- Colorless dipoles play 
two different roles in our approach. First, they are partons for the BFKL Pomeron. This role of 
the dipoles is not related to the large approximation. Instead of defining a probability to find 
Pomerons we search for a probability to find a definite number of dipoles at fixed rapidity. In this 
approach each vertex for splitting of one Pomeron into two Pomerons can be viewed as a decay of 
one dipole into two. Vise versa, merging of two Pomerons is an annihilation process of two dipoles 
into one. 

The second role of the color dipoles is that at high energies they are good degrees of freedom. 
This fact allows us to calculate splitting and merging vertices. However, we need to stress that the 
dipole model has been proven in the leading large Nc approximation only. The main assumption 
of this paper is that the dipole degrees of freedom can be in fact used for calculation of Pomeron 
vertices even beyond large Nc. 

Using this assumption we derive, in addition to previously known dipole splitting vertex r(l ^ 
2), two new vertices. The first one stands for the dipole recombination vertex r(2 —>■ 1). This 
vertex is derived by computing a lowest order loop diagram and it is essentially the same as the 
triple BFKL Pomeron vertex fT^ 120] . The only difference is in the normalization which allows us 
to use this vertex within the framework of the generating functional approach. 

The second new vertex is r(2 3), which accounts for the possibility of a dipole "swing". 
What we mean is that with some probability two quarks of a pair of dipoles can exchange their 
antiquarks to form another pair of dipoles. Naturally, this process has suppression. It is 

the vertex r(2 — 3) that correctly accounts for the Pomeron pairwise interaction in the BKP 
equation [23] and is absent in the usual form of the dipole evolution. 

We observe several advantages of our approach based on the generating functional. 

• Using the generating functional, we can separate the structure of the wavefunction of 
the produced dipole at high energy from rather complicated interaction of dipoles with the 
target at low energy. The latter are subject to non-perturbative QCD calculation and at the 
moment can be modeled only; 

• Our approach is based on colorless dipoles as the correct degrees of freedom at high energy. 
This fact allows us to have a natural description of hard processes in perturbative QCD, 
reproducing linear evolution equations, such as DGLAP ^26^ or/and BFKL [21] equations; 

• In our derivation of the linear functional equation we used a method which is closely 
related to the probabilistic interpretation of the Reggeon Calculus (see Refs. j221 123 I21j) ■ 
In doing so we establish a clear correspondence between the color dipole approach to high 
energy scattering, and Reggeon-like diagram technique providing a natural matching with 
high energy phenomenology of soft processes based on Pomeron. 

In the next section we describe the general formalism of the generating functional and its 
evolution equation taking into account both the emission of dipoles and their recombination. We 
also derive the equations for the scattering amplitudes which solve the problem of summation of 
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the BFKL Pomeron loops. In section 3 and in the appendix we discuss the dipole vertices for the 
process of transition of two dipoles into one. Pomeron interactions via two to three dipole decay 
is a subject of Section 4. Section 5 is devoted to study of dynamic correlations between dipoles. 
In conclusions we summarize and discuss our results. 

2 Dipole branching and Generating functional 

2.1 Classical branching process and equation for probabilities. 

We first consider a generic fast moving projectile whose wavefunction can be expended in a dipole 
basis. Note that contrary to many previous studies, we do not restrict ourselves to a single dipole 
as a projectile. 

= Y.an\n) (2.1) 

n 

Let us define a probability density Pn — to find n dipoles with coordinates ri, 6i, r2, 62, ■ ■ ■ ^'i, 
. . . r„, bji and rapidity Y in the projectile wavefunction. and bi are correspondingly the dipole's 
size and impact parameter, both are two dimensional vectors. We define P„ as a dimensionfull 
quantity which gives the probability to find a dipole with the size (from to + dvi). The 
integral 

J l[d\ Pn{Y -Yo;ri,bi,r2,b2,...rn,bn) = (2.2) 

and it gives the probability to find n-dipole with any sizes. This probability is conserved: J2n — 
1. 

Suppose that the following processes can occur as a result of one step in the evolution: 

1. The decay of a dipole with the size r and impact parameter b into two dipoles of the sizes 
r[ and ri^ and impact parameters b[ and , respectively: 

r(l ^ 2) = (r, b ^ r[, b[ + r'„ b',) ^ T,^2 (1; 1' 2') . (2.3) 

2. The annihilation of two dipoles with sizes ri and r2 and impact parameters bi and 62 into 
one dipole with the size r' and impact parameter b': 

r(2 1) = Ps^i (ri, bi + r2, 62 ^ r', b') = 12^1(1 2; 1'). (2.4) 

3. The process of interaction of two dipoles with sizes ri and r2 with a creation of one 
additional dipole (not factorisable to r(l — > 2) plus spectator): 

r(2^3) =r2^3(ri,6i + r2,62 ^ r[,b[ + r'^,b'^ + r',,b',) =T 2^3(1 2; 1'2'3'). (2.5) 
^We use notation r^, 6, for the initial state dipoles while r^, 6J for the final state ones. 
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4. The annihilation process of three dipoles into two dipoles: 

r(3^2) = r3^2(ri,&i + r2,62 +rs,h ^ r[,b[ + r'^,b'^) = Ts^2il 2 3 ; 1'2'). (2.6) 

In the remaining of this Section and in the next one we will focus on the first two vertices only. 
We will discuss 2 3 transition in Section 4. 



The equation for P„ obeys the classical branching process: 

i=l 



n „ 

-Ey dVfV^^2{i]l"2,')Pn{...nM ■■■) (2.7) 



dY 

+ E / dVi ri^2 (^'; t n) Pn-i (. . . r:, 6^, . . .) - E / dVf {t k ; 1') P„ (. . 



n-l 



+ T.fdV^ r,^i (z' k'; i) P„+i (. . . . . . r^, 6', . . .) 

ij^k'' 

Eq. ()2.7|) gives the general evolution for the probabilities P„. Eq. (j2.7|) must be supplemented 
by explicit expressions for the vertices r(l — > 2), r(2 1). By now, only the vertex r(l — > 2) 
has been calculated in Ref. pp. On one hand, all other vertices are formally suppressed by l/N^ 
and were omitted being considered as small. Moreover these vertices do not appear at all within 
the original formulation of the dipole model. On the other hand, it should be stressed that the 
Feynman diagrams which correspond to these vertices have the same order of magnitude as far as 
the Nc counting is concerned. For example, two diagrams in Fig. ^ show the Born approximation 
for r(l —>■ 2) (Fig. n]-b ) and for r(2 —>■ 1) (Fig. ^a ). They have the same suppression in N^. 
(^ but in the diagram of Fig. ^b this suppression could be absorbed in the amplitude of 

the interaction of two dipoles with the target while in Fig. ^a a factor has to be assigned 

to the vertex r(2 — > 1). 

Eq. ()2.7|) has a very simple structure. For every process of dipole splitting or merging there are 
two terms: the first one, with the negative sign, accounts for the probabihty P„ to decreases due to 
splitting or merging of one of n dipoles into dipoles of arbitrary sizes and impact parameters; the 
second term, with the positive sign, is responsible for an increase in probability to find n dipoles due 
to the very same processes. The first term includes the vertex r{n ^ k) integrated over the phase 
space of the final dipole: / dVj T(n ^ k) = J Yl'[ (Pr[ cPh'i V{n ^ k). For the second term we need 
to integrate over the phase space of initial dipoles: / dViV{n ^ k) = J Y[? (P'Ti(P'hiV{n k). 
Explicit expressions for the vertices F will be given in the next section. 

2.2 Generating functional and linear functional evolution 

The hierarchy ()2.7|) can be resolved by introducing a generating functional Z 

/n 
Pn{Y - Yo; ri,6i, r2,b2, . . . ,ri,bi, . . . ,rn,bn) Yl u{ri,bi) (f ncf k 

n=i i=l 

(2.8) 



5 



a) b) 

Figure 1: The lowest order diagrams for T {2 —>■ 1) (Fig.U\-a) and for r{l 2) (Fig.U\-b) . 

where u{ri, hi) = Ui is an arbitrary function of and h^. It follows immediately from ()2.7p that 
the functional ()2.8|) obeys the condition: at m = 1 

Z{Y - Fo; [u = l]) = 1. (2.9) 
The physical meaning of ()2.9|1 is that the sum over all probabilities is one. 

Multiplying Eq. ()2.7|) by the product HILi ^'^d integrating over all and 6j, we obtain the 
following linear equation for the generating functional: 

= x[u]Z[u]. (2.10) 

Let us introduce the dipole collective coordinate q with the integration measure q = d'^r d"^ b. 
The evolution kernel x is defined trough the operator vertices V 



X[u] = - J d'^qVi-,i{q] [u]) + j d"^ qi d"^ q2 Vi ^ 2 (qi, q2; [u]) 

Ai A2^2^2(gi;g2;N) + / AV2^i(g'; [u]) 
The functional form of the vertices are related to F's. 

Vi^i(g,M) = I dUid\2Ti^2{q ^ qi + q2)u{q)j^^. (2.11) 

Vi^2{q[: q'2: N) = / d\V,^2 {q^ q[ + q'2) <q'i) ^(92) ^ • (2-12) 

V2^2 (gi; g2; M) = \ I d^qT2^iiqi + g2 ^ g) u{qi)u{q2) ^ ^ . (2.13) 

2 J ou[qi) du{q2) 

V2^i{q;[u]) = \ I dUidU2T2^i{qi + q2 ^ q) ^^l) ' (^.14) 

2 J ou[qi) du{q2) 



The functional derivative with respect to u{q) = u{r, h) plays a role of an annihilation operator 
for a dipole of the size r, at the impact parameter h. The multiplication by u(r, h) corresponds to 
a creation operator for this dipole. 

Eq. ()2.10|) exhibits a quantum mechanical-like structure with the operator x viewed as a 
"Hamiltonian" of the evolution. In Eq. ()2.1H1 x is constructed in terms of dipole creation and 
annihilation operators. The evolution operator x describes a two dimensional fully quantum 
(nonlocal) field theory of interacting dipoles. 

2.3 Evolution of dipole densities 

The n-dipole densities in the projectile Pniji, bi, . . . , r„, 6„) are defined as 

1 6 

Pnin, bi ... ,rn,bn) = -^Y[t—Z{[u]) \u=i (2.15) 

1=1 ' 

Differentiating Eq. ()2.10j) n times with respect to u we can obtain a hierarchy of equations for p„ 
Refs. ini 12 E] and rewrite Eq. ()2.7p in the form: 



dpn {ri,bi] . . .■,rn,bn]Y) 
dY 



n— 1 „ 

^ / rf^gri_2(g; qiqn)pn-i (. . .g. . .) (2.16) 

i=i 

n „ n „ 

+ 2^ d^q' d^q Ti^2{,q'] qqi) Pn{---q' ■■■) - E / ^^^i ^^^2 ^i^2{qi] qW2) pn{,- ■ ■ qi ■ ■ ■) 
i=i •' 1=1 •' 

n „ 

- H j d'^qT2^iiqtqk; q) Pn{. ■ ■ ,qi,- ■ -qk ■ ■ ■) 

i,k,ij^k 

n . n . 2 

- 2 ^ J d'^q d^gTs^i [qqi] q') p„+i {. . . qi . . . q) + ^ j \{ d^q[ Ta^i [q'^q'^] qi) pn+i {q[ ■ ■ ■ q'2) 



i=i •' i=\ 



Eq. (j2.16|) presents a general structure with so far arbitrary vertices r(l — 2) and r(2 — 1). 
Explicit expressions for F will be presented in the next section. The diagonal part of the evolution 
due to the vertex r(l 2) is the large Nc limit of the BKP equation ^25J in coordinate space, 
n = 1 corresponds to the evolution of the BFKL Pomeron. 

The vertex r(2 1) (Eq. (j3.3(Jj) ) obeys the following property / d^qV2^ \ {qiqk, q) = ^. 
Consequently, it is only the last term proportional to r(2 ^ 1) who survives in Eq. (j2.16p . 

2.4 Scattering amplitude 

As was shown in Refs. 13] , the scattering amplitude is defined as a functional 

N {Y- [7,,]) = -Y: / &i; . . . r„, Fq) 11 (F, [u,]) \u^=id\ d\ 

n=l i=l "'^i 



^We are thankful to our referee who actuaUy noticed this property. 
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oo „ n 

= -E(-l)" / Inin, bi;...r^,bn;Yo)pn in, bi;...r^,bn;Yo) II d\d\. (2.17) 

n=l 1=1 

The amplitude for simultaneous scattering of n dipoles off the target is denoted by 7„. It has to be 
specified at the lowest rapidity (Yq). Using the ansatz {(Y, [7^]) = (71(F), 72(1^) . . . 7i(F) . . .) 
and having 



^ = y / rr d'r,d%, — (2.18) 

dY J l\ ' ' 6jn dY ^ ^ 



we can recast the hierarchy of equations (j2.16j) into the Balitsky-type chain for the scattering 
amplitudes (see Ref. 0) 

dY - ^^-^^^ 

= 2^ / d^q' d'^qri^2{qi; qq')in{- ■ -q' ■ ■ ■) - I d'^q'id'^q'2^i^2{qu q'iq'2)in{- ■ ■ ,qi- ■ ■) 

i=i i=\ 

n-l 



- H j d'^qd^q'Ti^2{qh qq') In+i {■ ■ ■ q ■ ■ ■ q') + J d^ qT2^i {qiqj] q) -fn-i {qi ■ ■ ■ qj ■ ■ ■ q) 

i=l i^j 

n . n . 

- 2 y / rf^g d'^q'r2^i [qqu q') In-i {...qi...q)+ Yl / ^2^1 (^^^^5 g) 7„ (. . . . . . . . 



i=l i,k,i=^k 

The great advantage of Eq. ()2.17p is the fact that this equation allows us to take into con- 
sideration in the most economic way the interaction of low energy dipoles with the target. For 
example, assuming pi(lo) = 5*^^^ (r — ri) 5^^^ (b — bi^ while p„>i = 0, we obtain that the total 
amplitude of a single dipole scattering equals to 

N{r,b;Y) = ^i{r,b;Y) . (2.20) 

If we assume the projectile be built out of two dipoles with 

p2(Fo) = S^^^ {r - ri) 6^^^ (b - 5(2) (f" _ fa) 5(2) (b' - 62) and pi = 0, p„>2 = then 

N {r,b;r',b';Y) = 72 (r, 6; r'6'; F) . (2.21) 

Eq. (j2.19j) is an evolution hierarchy for dipole amplitudes. Apparently it involves loop processes. 
The equation is most general for postulated vertices. It has a very similar structure as suggested 
by lancu and Triantafyllopoulos in Ref. (Eqs.(6.6) and (6.7) of this paper). In the following 
section we will present explicit expressions for all F. The exact vertex F(2 — 1) found by us 
(see Eq. (|3.3U|) ) does not seem to coincide with the vertex suggested in Ref. ^3 ^^J kinematic 
region^ . 

In fact, the last two terms in Eq. (j2.19|) vanish for the vertex F(2 1) given by Eq. (j3.3(jp . 
Nevertheless we prefer to keep these terms explicitly in the hierarchy Eq. ()2.19p . The only reason 

•^Our vertex does coincide with the one found by the authors of in their paper |2H1, which appeared after 
our preprint started to circulate. 



8 



behind keeping them is that in practical apphcations one may attempt to approximate or simphfy 
the vertex r(2 —>■ 1). For an approximate vertex the last two terms might not vanish. Diagram- 
matically these terms are part of the 4 — > 4 reggeized gluon transition which may contribute for 
some kinematics where the underlying probability conservation is important. 



3 Dipole vertices 

3.1 r(i 2) 

The vertex for the decay of one dipole into two has been derived in Ref. pQ: 

T,^2{r,b r[,b', + r',,b',) = (3.22) 
^s:£^ 5^'' (r - K - K) 5^'' {b[ - b + Irl) 5^^^ (b', - b - K^) . 

As has been discussed, this vertex leads to reduction of the probability to find n dipoles due to 
decay into two dipoles of arbitrary sizes. This reduction is related to 



r.2 

ip r'^ [r — r' 

where as = asNc/ir and p is the infrared cutoff. The growth term is proportional to 



— 2 

«5^(r) = / dV'T^^2{r,b r[,b', + r',X) = / ^J^rf^^ r' (3.23) 



J dVT,^2ir,b ^ r[,b'^ + r'^X) = J d^r d^b T,_,2 ir,b ^ r[,b'2 + r'^X) (3-24) 

5^'^ {h\ - A - \{r, + r1)) . 



2 TT rf r'i 



So far, most of the discussions of dipole evolutions were bounded to the vertex r(l — > 2) as 
all other vertices are l/N^ suppressed and were considered as small corrections. Below we will 
include several new vertices, which are of the order These additional vertices give important 

contributions in the deep saturation region (see Refs. jHl 1121112111511121 for more detailed discussion 
of this subject) and hence need to be accounted for. Though we do not pretend to be able to 
accommodate all of the A^^^ corrections using dipole degrees of freedom, we believe the contributions 
we aim to include are dominant at high energies. 



3.2 r(2 1) 

To find the vertex r(2 —>■ 1) we analyze the first enhanced diagram shown in Fig. |21-a. As was 
shown in Refs. ^21 ^] the expression for this diagram is 

J d^ ri d^ bi d^ r2 c/^ 62 d^ r[ d^ b[ d^ r'^ d^ b'^ P2 {r[ b[; r'^,b'^) 

P2{ri,bi;r2,b2) -f^s^ {ri,r[,bi - b[^ Tba (^2, ^2, &2 - ?) (3.25) 
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where 7^]^ is a dipole-dipole elastic scattering amplitude due to exchange of two gluons. The 
expression for this amplitude is well known (see Refs. [2^ I19j ) 



7« fn r'-l - b') - In^ (A& + R? {Ab - Rf , ^ 



a 



S i„2 



In" 



2 ^ ^xi-x[)'^{yi-y[]'^ 



32iV2 V(^i-^i')' (^i-^i)V 
where S = |(ri + r^), R = |(ri — r']'), Ab = bi — b[] and ri = Xi — yi, r[ = x[~ y[. 

Eq. (j3.25|) follows from the fact that we can view this diagram in the following way. The upper 
dipole (xio in Fig. |2)) evolves with normal vertex r(l 2) in the rapidity interval y — yi (see 
Fig.|21-a) while the low dipole (xi'o' in Fig. El) also evolves but in rapidity interval y2. As the result 
of these two evolutions there are two dipoles with the sizes ri = 0:20 and r2 = X12 at the rapidity 
yi and two dipoles with sizes r[ = X2'o' and r'2 = X112' at the rapidity y2- Each pair elastically 
rescatters by the exchange of two gluons leading to Eq. ()3.25j) . 

Alternatively, Eq. ()2.7p gives another expression for the same diagram 

I (f n (f bi (f r2 (f b2 (fr (fbP2 {r[ b[ r'2 b'2) Fs^i (ri, 61 + ra, 62 ^ r, b) 7^ (r, r'; b - b') 

(3.27) 

where r' = xi'o' in Fig. |21 Comparing Eq. ()3.25|) and Eq. ()3.27|) we obtain the following equation 
for r2^i (ri, 61 + r2, 62 — * r, b): 

j SrdPbV 2^1 {nM + r2,b2^r,b) -i^B A {r.r\b -b') = (3.28) 

r[ b[ d' r', d' b'2 T,^2 (r', b' ^ r'„ b', + r'„ b',) 7^ {n, h - b[) 7^ (r2, r'2, h - fc^) . 

Eq. ()3.28|) is the basic equation from which the vertex r(2 1) can be extracted. To this goal we 
need to invert Eq. ()3.28p by acting on both sides of it by an operator inverse to 7^^ in operator 
sense. Fortunately, this operator is known to be a product of two Laplacians: 

A.A,7S(x,l/;a:',y') = al [5^'^ [x - x') 5^'^ [y - y) + 5^'^ [x - y) 5'^'^ (y - x')) (3.29) 

with X = b + r/2 and y = b — r/2 being the coordinates of quark and antiquark in the dipole (r, b). 
For the vertex r(2 — 1) we finally obtain 

r2^i (ri, 61 + r2, 62 ^ r, 6) = ^ A, Ay f d^ r[ d^ b'^ d^ r^, d^ b'^ Ti^2 (r, b ^ r[, b'^ + r^, b'2) x 

as J 

X 7BA(n,ri,6i - b[) lBA{r2,r'2;b2 - b'^ . (3.30) 

In Appendix A we present a method for evaluation of the expression ()3.30|) . We arrive at the 
result given by Eq. Since the general expression is rather complicated we consider now 

simplified estimates valid in several different kinematic regions. 
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In the region where Ab <^ r' and r <^ r' Eq. (|3.26|) leads to a simple expression 

7a(r,r';A6) = (3.31) 

Using Eq. ()H.H1|1 we can rewrite Eq. ()H.28j) in the simple form if we are looking for the contribution 
in the following kinematic region: 

ri r2 r'^ ^ r'2 ^ r' and r ^ r' 

In this kinematic region using Eq. ()3.22|) and Eq. ()2.7|) we obtain that the r.h.s. of Eq. ()3.28|) is 
equal to 

r.h.s. of Eq. (jS2HD = fej J d^r[d^ r'^d^ b[ d\ (3.32) 



^(^) (f ' - K - K) s^^^ {b[ - b + If;) ^(^) (fei - b - Ik) 



l,iV2 j y r;2r;2 ^2^2 «5 ^^2 y ^2 

The l.h.s of Eq. fl3.28|) has a form 

l.h.s. of Eq. (IT^ = (3.33) 

^2 J d'^rd'^br2^i{ri,bi+r2,b2^r,b)(^BA{r,r';b-b"j 
32 a\ 



TT r' dr T2^i{ri,bi + r2,62 ^ r,b) 



In the last equation we used the fact that in Eq. ()3.2(ip the typical b is about the size of the large 
dipole {b ^ r'). Substituting Eq. (HO^ and Eq. S^H^i in Eq. (HT^ we obtain T (2 ^ 1) in the 
form 

(4«5f 1 



r2^i(ri,6i + r2,62 ^ r, 6) = - (3.34) 



Repeating the same calculation but in a different kinematic region: 

ri ~ r2 ^ r\ ^ r'^ and r'^ ~ ^"2 ^ while r <^ r' 

we obtain 

r2^i(ri,fei + r2,62 - r,6) = ^ (3.35) 

The full expression for F (2 1) is rather complicated as can be seen from Eq. (IX22ll . In 
corresponding limits it reproduces Eq. ()3.34|) and Eq. ()3.35p . 
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a) b) c) 

Figure 2: The first enhanced diagram and the vertex of annihilation of two dipoles into one dipole. 

4 Pomeron interaction: 2 — > 3 transition vertex 

In this section we further extend the dipole model by introducing an additional splitting vertex, 
r(2 — >■ 3). Our main goal here is to account for Pomeron pairwise interactions via exchange of 
a single gluon. For the first time this process was included in the double log approximation of 
pQCD in Refs. [HHIIS]- It is also most naturally included in the BKP equation providing 
corrections to the dipole evolution discussed above. 

The inclusion of the above Pomeron interactions in terms of dipole degrees of freedom is not 
a straightforward task. We face two problems here. First, the contribution we are looking for is 
a process in which a gluon is emitted (in the amplitude) by one dipole and then reabsorbed (in 
the conjugate amplitude) by another dipole. This is an interference contribution, which does not 
admit a probabilistic interpretation. The exact expression for Nc corrections to the dipole evolution 
known^ from Refs. [321 iSSj can, nevertheless, be projected onto dipole degrees of freedom. By 
introducing the vertex r(2 — 3) we take into account only the diagonal contributions factorisable 
in terms of dipoles. We trust that the rest of the Nc corrections contribute to multi gluon t- 
channel states only. The latter, n-gluon states are known to have intercepts smaller than that of 
n/2 Pomerons and thus could be ignored at high energies. 

The second problem is in the fact that dipoles are natural degrees of freedom in the large 
Nc limit only. Beyond large N^, the dipole basis ()2.H) is overcomplete. In particular, a single 
space configuration of two pairs of quarks and antiquarks can be counted twice as two different 
pairs of dipoles (provided all quarks are mutually in a color singlet state). As a result of working 
with overcomplete basis there will be a nontrivial overlap between probabilities to find a different 
number of dipoles. 

Having sorted the above problems out, we propose the following vertex x^^^ to be added to 
''We thank Yu. Kovchegov who drew our attention to Ref. Q2i after our preprint started to circulate. 
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the dipole evolution kernel x- 



= / T2^3ix2,Xs,X,) X (4.36) 

Jxi,X2,X3 X4, X5 

X X 

X [1 - m(xi,X4)] [u{x2,Xs) - U{X2, X5) U{X5, X3)] 



6u{qi) 6u{q2) ' 

The operator ()4.3(i|l describes the following process. First, it annihilates two dipoles qi = (xi, X2) 
and q2 = {xs, X4) (Fig. El). Then the dipoles are regrouped into 14 (spectator) and 23. The latter 
subsequently decays through the usual 1 — 2 dipole splitting process {u{xi,X4^) {u{x2,X3) 
— m(xi, X4) ■u(x2, xs) ^(xs, X3)) term in the operator). The vertex r(2 ^ 3) has the usual dipole 
splitting form (r(l 2)) suppressed by a factor A*"^: 

12^3(2;, z) - 



2 7riV2 {x - zf{y - zf 

In fl4.3fi|l we have also subtracted a term with the spectator u{xi,Xi) set to unity. This 
subtraction is needed to remove the double counting: the decay of a single dipole (23) has been 
already accounted for in the normal 1^2 dipole evolution. This subtraction can be also thought 
of as originating from the overcomplete basis we are dealing with. We will find below that this 
subtraction is crucial to prevent the operator (j4.36|) from generating Pomeron loops at the level 
of scattering amplitudes ^. 

For the evolution of the dipole densities pn the operator x^^^ generates the following contri- 
bution (for n > 2) 

gy = 1^ T2^3{Xj,yk,yj) Pn~-i{- ■ ■Xi,yk- ■ ■Xj,yi. . .)d {Xk - Vj) 

i,j,k=l,i^j^k 
n „ 

+ J ^-i^sixj, y, yj) pn{- ..Xi,y... Xj, y,...) 

n „ 

j (fz r2^3(xj, Vj, z) pn{. . . Xi, yj... Xj, Vi . . .) (4.37) 



The evolution of the scattering amplitudes receives additional terms (for n > 2 
din (gi; . . . ;g„;F) 



BY . . ... 



J d^z r2^3{xj,yi, z) [7„(. . . Xi, yj... Xj, z . . .) + (4.38) 



7„(. . . Xi, yj . . . z,yi. . .) - 7„(. . . Xi, yj . . . Xj, yi . . .) - 7„+i(. . . x,, yj . . . Xj, z . . . z,yi . . .)] 

Eq. ()4.38p supplemented by the usual dipole evolution generated by the vertex r(l — > 2) is 
believed to be a very good approximation to the Balitsky-JIMWLK evolution. The advantage 



'^We have missed this subtraction in the first preprint version of this paper. We are most thankful to our 
coheagues Ian Balitsky, Jochen Bartels, Al Mueller, Yura Kovchegov, Alex Kovner, and the referee whose criticism 
helped us to solve the problem. 
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b -1/2( X2 -X4) 



b + 1/2( x^-X3) 



^5 



X^2(>'i) + X34(r2) 



I *-X25(r'2) + X35(ri) 



Figure 3: The process Xu + X34 + X35 + X45. 



of our formulation is that it is given entirely in terms of dipole degrees of freedom. We will 
demonstrate in the following section that the above evolution happens to coincide with the one 
found in Ref. pij by analyzing corrections arising from the QCD triple Pomeron vertex ^3]. 

Finally let us comment about 3^2 transition vertex. The process of 3 — > 2 has 
compared to the leading 1 — > 2. Indeed, it progresses in two stages: the first one is the annihilation 
of two dipoles into one. Such a process has suppression. Then two remaining dipoles "swing" 
quarks (see Fig. O)) and this has an additional l/N^ suppression. Therefore, r(3 2) is of the 
order of l/N^ and will be neglected. 

5 Nc correlations due to 2 — > 3 vertex 

Let us combine Eq. (j2.19|) and Eq. (j4.38p but neglect the vertex r(2 1). We would like to 
find a procedure which would allow the equations entering the hierarchy to decouple from each 
other. In case of original Balitsky's hierarchy this was achieved by assuming absence of target 
correlations which means substitution of the Kovchegov's factorization [H]: 

n 

7„(ri,6i r„,6„;F) = J| 7^ (r^, 6^; F) . (5.39) 

i=l 

The whole hierarchy respected the factorization leaving only one single equation (BK) unresolved. 

Since we have dynamical correlations, the hierarchy of Eq. ()2.19p +Eq. ()4.38|) obviously does 
not admit the factorization of Eq. ()5.39p . An intuitive solution would be to introduce pairwise 
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correlations which would hopefully reduce the hierarchy to two coupled equations. The natural 
generalization of Eq. ()5.39|) is to introduce two-dipole correlation C iji, bi;rk,bk; Y) in the form 

7„(ri,6i r„,6„;r) = (5.40) 

n n n 

n 7i(r.,6i;r) + n C{r,,k;rk,h;Y) ^,in,bi;Y) 

1=1 i=i,fc=i,j^fc/=i,/^ij^fc 

Eq. ()5.4()|1 can be written compactly by introducing the operator Q, such that 

n 

-fn{ri,bi . . . , rn,bn]Y) = Ql'Ji] Yl^iiri, bi) , (5.41) 

i 

with 



g[7i] = Exp(^J d%d%Ciq^, q^; Y) 



5 5 



57i(gi) 57i('?2)/ 

We have checked that, though the introduction of correlations in the form Eq. ()5.40p is very plau- 
sible idea, this ansatz does not make the system of hierarchy equations to decouple. Nevertheless, 
we can try to estimate the influence of the new vertex by taking into account the correlations 
between dipoles in perturbative way considering them small. To this goal we will focus on the 
first two equations of Eq. (|2.19p +Eq. (j4.38|) which will allow us to determine the evolution law for 
7i and C. Introducing K as the usual dipole kernel 



K{x,y;z) = (5.42) 



the equations for 71 2 read 



dY 

di2 ixi,yi;x2,y2) 



as / K {x, y; z) (- 71 (x, y) + 71 (x, z) + 71 (y, z) - 72 (x, z; y, z)) 



(5.43) 



dY 
2 

as 



Y K {xi, yi] z) [72 (xi, z- Xj, yj) + 72 (yi, z; Xj, %•) - 72 (x^, yf, Xj, yj)] 
•^^ i,j=i,i^j 

+ ^ / ^i^v Vi'^ ^) Xj, z) + 72(xi, yj\ z, yi) - 72(xi, yj\ Xj.yi)] 

In Eq. ()5.43|) we have omitted terms proportional to 73. Substituting 

7i(x,i/;F) = N{x,y-Y) 

and 

72(xi,?/i;x2,i/2;>^) = A^(xi,?/i;F) N {x2,y2\Y) + C (xi, X2, 1/2; >^) 
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we obtain assuming the correlation function C is small, C < N"^: 

dN{x,^,Y) ^ as J d'^z K {x,y;z) {Af {x,y;z) - N {x,z) N {x,z) - C {x,z;y,z;Y)) 



= as 

where we define 



(fzK{x,y;z) (Af{x,y;z) - C {x, z;y, z;Y)) ; (5.44) 



N'{x,y;z) = N{x,z) + N{y,z) - N {x,y) (5.45) 
Af{x,y;z) = Af{x,y;z) — N{x,z) N{x,z) 
C{xi,yi,z;x2,y2;Y) = C {xi, z;x2,y2;Y) + C {yi, z;x2,y2;Y) , - C {xi,yi;x2,y2;Y) 
The equation for 72 becomes an equation for C 

dC {xiyi;x2,y2;Y) 



dY 



(5.46) 



+ 



+ 



(/ . (K C ... Y)^ K (.„ C x„ Y)) 

d"^ z K {yi,X2]z) {N{xi,yi) + N {x2,y2) - N{xi,y2)) {j\f {yi,X2-, z) + M {xi,y2; zfj 



as 
2 

as 
2 m 



It is interesting to notice that Eq. ()5.4fi|l can be reduced to the equation of Bartels, Lipatov 
and Vacca Indeed, we can introduce a new function AC: 

AC{xi,yi-X2,y2]Y) = (5.47) 

C(xi,yi;x2,i/2;y) - ^ {N{xi,yi-Y) + N {x2.y2;Y) - N {xi,y2]Y)f 

Using Eq. ()5.44|) we can reduce the set of Eq. ()5.44|) and Eq. ()5.46|) to a different set of equations, 
namely, 

dN{x,y-Y) 



dY 



(5.48) 



as J d'^ z K {x,y;z) (^N'{x,y;z) - N {x, z) N {x, z) - ^^Af^ {x,y; z) + AC {x, z;y, z;Y)^ ; 

dAC{x,z-y,z;Y) 



dY 

as_ 
2 



(5.49) 



+ ^ I <Pz {K{yi,X2:z)N{yi.z:Y) N(x2.z^,Y) + K (xi,y,; z) N (y„ z^,Y) iV(i,,2;F)) 



C 
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These two equations are the same as were proposed by Bartels, Lipatov and Vacca j^lj. Our 
derivation suggests also a physical meaning of the modified Balitsky-Kovchegov equation (see 
Eq. ()5.48|) ). The Balitsky-Kovchegov equation is a mean field approximation while Eq. ()5.48p 
takes into account the correlation related to possibility for grouping of two dipoles in a different 
way with suppressed probability. Therefore, it plays a role of Fock term in Hartree-Fock approach, 
which is a natural next step in the mean field approach. AC is a real dynamic correlations which 
as one can see from Eq. ()5.49|) grows with energy. We have neglected terms of the order ACN in 
comparison with N'^ -term. Therefore, we can trust the equations only for AC < N. For higher 
energies we need to develop a more general approach. 

6 Conclusions and Discussion 



In this paper we have extended our linear operator approach applied to dipole evolution. The 
evolution kernel x can be viewed as a "Hamiltonian" of the evolution. It is constructed in terms 
of dipole creation ad annihilation operators. By introducing the recombination vertex F (2 — > 1), 
the evolution operator has been promoted to a fully quantum two dimensional field theory of 
interacting dipoles (Pomerons). 

The main resuhs of this paper are Eq. (jSHDl), Eq. (jTTT^ . Eq. ^U^ . Eq. and Eq. 

supplemented by the explicit expressions for the vertices F (2 — 1) and F (2 — 3) which both 
are proportional to the second functional derivative with respect to Ui. Our approach is an 
extension beyond the Balitsky one j3] , based on the Wilson loops, as well as beyond the Color Glass 
Condensate approach (JIMWLK equation ^01 E])- Though the JIMWLK equation takes into 
account all corrections, and which are only partially accounted for by the vertex F(2 —>■ 3), 

they do not include the recombination vertex F(2 — > 1) which is a major step beyond this equation. 

We have accounted for dynamical correlations that stem from possibility of merging of two 
BFKL Pomerons. It is illustrative to consider a simple toy model in which we assume that 
interactions do not depend on the dipole sizes (see Refs. Pfl |21 CS] for details). The master 
functional equation (see Eq. ()2.10|) ) for this model degenerates into ordinary equation in partial 
derivatives 

W = -^^^ ^ 2) nil -u)^ + Ti2^ 1) ^1 - ^ + r(2 - 3).. (1 - uf^^ (6.50) 

We can introduce a generating function for the scattering amplitude using the relation 

N{Y;u) = 1 - Z(r;l-7) (6.51) 

To obtain the scattering amplitude we need to replace 7 in Eq. ()6.51|) . by the amplitude of 
interaction of a dipole with the target. For Eq. ()6.50|) can be rewritten in the form: 

B N f) N f)'^ N r)"^ N 

^ = F(l ^ 2) 7(1 - 7)^ + r(2 ^ 1) 7(1 - 7) ^ + r(2 ^ 3)7^ " ^ W ^''''^ 

if 7 is small we can reduce Eq. (I6.52|l to a simpler equation 

dN ON 

-^ni^2),-. (6.53) 
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The solution of Eq. (j6.53|) is a Pomeron with the intercept r(l 2): 

The rest of the terms in Eq. (|6.52j) are responsible for Pomeron interactions (see Fig. 0}. 

P P p 




Figure 4: Pomeron interactions described by Eq. 

As we see from Fig. E] two vertices r(2 — > 1) and r(2 — > 3) are responsible for different 
processes of Pomeron interaction. At first sight r(2 — > 1) is much smaller than r(2 3) and 
can be neglected. However, we can make such a conclusions only if we will find out what value of 
u { or j) are essential for high energies. Therefore, the vertex r(2 — >■ 1) can be still relevant in 
certain kinematic domains. To answer this question we need a detailed analysis of Eq. (|2.19|) and 
Eq. (j4.38p which is beyond the scope of this paper (this question is addressed in Ref. jSS])- 

Fig. presents some examples of Pomeron diagrams which correspond to different approaches 
that has been discussed in the past: the GLR equation |H] (see Fig. |Sl-a) which, in our approxima- 
tion, coincides with the BK |SJ [H] and JIMWLK fTT] equations; the lancu-Mueller approach J2] 
(see Fig. |S]-b). Fig. |SJ-c shows a typical diagram that can be incorporated using Eq. ()6.52|) . Finally, 
in Fig. El-d we plot the diagrams that one needs to sum in order to reliably consider nucleus - 
nucleus interactions. In general such diagrams are difficult to sum, but we have an experience 
that in the simple model of Eq. (j6.52|) . this summation can be performed jSUSE]- 

It is important to stress that by introducing the vertex r(2 — 3), we have taken into account 
only the leading Nc corrections. For n dipole densities with > 2 we should have color correlations 
which cannot be presented in the dipole basis. We believe, however, that these correlations are of 
no significance at high energies. 

The importance of correlations have been already noticed in Refs.jSOl 1^ Eq- ()6.52j) 

illustrates a complexity of the problem since even this oversimplified equation has not been solved. 
The expansion in correlations allows to shrink the infinite hierarchy of equations to a system of 
two coupled equations. This reduction provides a method for estimating importance of both the 
Nc correlations and Pomeron loops. We demonstrated that correlations should be essential at 
high energies and suggested a consistent approach to take them into account. 
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GLR(JIMWLK) lancu - Mueller 




Figure 5: The typical Pomeron diagrams for interactions described by Eq. \6.5^) . Fig.\^a describes 
the GLR approach /5/ which for Eq. \6.5i^) coincide with the Balitsky-JIMWLK approach ni[ 
Fig. \^b corresponds to lancu-Mueller approach]^ which suggests the way out of the JIMWLK 
approach and can be justified in limited region of energy. Fig. \^c shows the general type of the 
diagrams that can be summed in the framework of the approach based on Eq. \6. 5^) . Fig. \E^d are 
diagrams that we need to sum for nucleus-nucleus interaction at high energy. 
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In this paper we have considered a merging process of two Pomerons into one only. In general, 
there exist higher order processes accounting for a possibility of many Pomerons merging into one. 
A formal resummation of these processes has been reported in recent Ref. [SZ] and also in Ref. 

m- 

We hope that we propose the simplest way of dealing with the Pomeron loops which is equiv- 
alent to the reggeon calculus for BFKL Pomeron but has an advantage of clear probabilistic 
interpretation in the rest frame of one of the colliding particles. We hope that clarification of all 
assumptions in our approach will lead us to deeper and more transparent understanding of physics 
in the saturation domain. 

Finally, let us comment on two recent papers jHHlEHl which appeared practically simultaneously 
with ours and contain features close to those presented here. In fact the formal expression for the 
vertex r(2 — > 1) (Eq. is identical to the ones of Refs. [HSIEH]- This equivalence has been 

proven in a later Ref. j^lj. The main difference is that we have extended the method. Apart from 
giving the formal expression for the vertex r(2 1), we also introduced a formalism needed for 
its evaluation (see Appendix). At the end, we were able to obtain a first analytical evaluation 
of the vertex bringing it to the level ready for computer simulations (Eq. ()A.22j) ). The diagonal 
transition 2^2, which guarantees probability conservation, vanishes if the exact expression for 
the vertex r(2 —>■ 1) is used. Generally this term should not be neglected if an approximate 
vertex is used for practical applications. In addition we have included the vertex r(2 — > 3) in our 
consideration. 
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Appendix A: Calculation of F (2 — > 1). 

In this appendix we find the solution to Eq. ()3.28|) . Our approach is based on the main properties 
of the BFKL kernel which have been studies in details in Refs. P^ I2Uj. First, we rewrite Eq. ()3.26p 
in the form of the contour integral over h ^3 120] , namely, 

1^"^ {x,y]x',y') = ^{^BA{x,y;x',y') - ^BA{x,y;y',x')} (A.l) 

where 

^BA{x,y;x',y') = (A.2) 
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■'jBA{h;x,y;x ,y) - 



h 



Ja-ioo 27rih^ 16 Ja-ioo 27ci \{x — x'Y {y — y'Y ) 

where x, x' and y' is the coordinate of quarks and antiquarks in the interacting dipoles with the 
size r = X — y and r' = x' — y' . Introducing complex numbers instead of vectors x = Xi + ixi and 
X* = X\— 1X2 {xx* = (xY), we can rewrite ^ba {h; x, y; x', y') in the form: 



(A.3) 



At first sight Eq. (jA.lll does not lead to Eq. ()H.2(j|l . Indeed it gives 

^^^{x,y-x\y') = (A.4) 



32 \{x — x'Y {y — y'Y ) \{y — x'Y {x — y'Y {y — y'Y [x — x'Y ) 

The replacement of the Born amplitude Eq. (j3.26p by Eq. ()A.1|) is a major step for what follows 
and has to be justified. We refer here to the work of Lipatov 19 who showed that the Born 
amplitude could be written in the form of Eq. ()A.1|) (see Eq.llO of the first paper in Ref. |19j). 
The main idea of Ref. jTHj is that two expressions Eq. ()A.1|) and Eq. ()3.26|) lead to the very 
same results if used for calculations of physical observables (for example 7 — 7 scattering). Both 
expressions satisfy Eq. (j3.29p and hence they differ by a function ^, which does not depend on one 
of the coordinates x (or y). Lipatov showed that, thanks to the properties of the impact factor 
$(x, q) (see Eq. 109 in Ref. [IH]), the integral over x (or y) of the impact factor convoluted with 
^ vanishes. This property of the impact factor implies that a function, which does not depend on 
one of the coordinates, gives zero contribution to any physical process. Moreover, the well known 
BFKL Green function j2I] was calculated using Eq. ()A.1|) as initial condition. To be consistent 
with the use of the BFKL kernel, Eq. flA.l|) has to be taken as the Born approximation. We will 
see below that this replacement allows us to evaluate the vertex r(2 ^ 1) (Eq. ()3.30|) ). 

In what follows we deal with the first term in Eq. ()A.1|) but it is a trivial algebraic exercise to 
obtain a result for the full Born amplitude of Eq. ()A.1|) . 

The r.h.s. of Eq. ()3.28|) we rewrite, using Eq. ()A.2|) and Eq. ()A.3|1 in the form 

ra+ioo dhi 1 ra+ioc d h2 I f , , ^ f X - y \ ( X - y 

as / : -pr / : -pr I dz dz 



a-ioo 47ii hi Ja-ioo Aixih^J \{x — z) {jj — z) ) \{x — z) {y — z) ) 

^BA{hi;x,z;xi,yi) ^BA{h2; z,y;x2,y2) (A. 5) 

The integrals over z and z* can be computed using formula 3.211, 9.182(1) and 9.183(1) of 
Ref. jini- Indeed, Eq. ()A.5|1 can be rewritten as follows 

rv?, ra+ioo rlh, 1 ra+ioo rl h „ 1 

«5 [ a til I f ati2 I ^^ g^ 



16 A^^)^ Ja-ioo 4:7Ti J a-ioo 4 71 i hi 
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Using the notation x — y = r and x — z = r', we have 

r / r'{xi-yi) \ / jr' - r) (xa - ^2) \ 
(r-r')r' \{x - Xi) (x - r' - yi) J \{y - yi) {x - r' - X2) ) 

We calculate J dividing the integration over r' in three regions^, namely, 

I = h (-00 < r' < 0) + /a (0 < r' < r) + J3 (r < r' < 00) (A.8) 
Let us first compute the integral from the second region 



dr' 



r'{xi-yi) f {r' - r) {x2 - y2) 



(r — r') r' \{x — Xi) {x — r' — yi) I \iy — 2/2) {x — r' — X2) 



( {xi -yi) {x-y) \^^ f {x2 - y2) [x - y) \^'' 
\{x - xi) {x - yi) ) \{y - y2) {x - X2) ) 



(A.9) 



(A.IO) 



T{hi + h2) V ' ' ' ' {x -yiV {x -X2) ) 

( {xi - yi){x - y) V'^ ( {x2 - y2){x - y) V'^ 
\{x - xi) {y - yi) ) \{x - X2) {y - y2) ) 

v{h,)v{h2) p. (, , , i^-y) ix2-yi) \ 

" Tih, + h2) ^FAh,,h2M + h2, (,_,^)(^_^^)J 

_^ {hi + h2) f {xi-yi){x-y) \^' f {x2 - y2){x - y) \^' (All) 
hih2 \{x - xi) {y - yi) J \{x - X2) {y - y2) ) 

Eq. was obtained using 3.211 of Ref. [10] • Obtaining Eq. (jA.llI) we take into account that 

only small values of hi and /12 will contribute to the integral of Eq. ()A.5j) . 

In Eq. ()A.9j) and Eq. (jA.llj) Fi and 2F1 = F denote the hypergeometric functions (see formula 
9.10 and 9.180(1) in Ref. [40j). 

J3 is the integral of Eq. ()A.7|) for r' > r, namely, 



r'{xi-yi) f (r - r') (X2 - 1/2) 



r (r — r') r' \{x — Xi) {x — r' — yi) ) \{y — y2) {x — r' — X2) ) 

hi / f ^ \\h2 

(A.12) 



\xi -yi) \ ' ( {x2 - 1/2) \ 
\x-xi) ) \{y- 2/2) / 



r(l)r(M / {y - yi) {y-y2r 

V{l + h2) \ [x-y) [x-y) 

_^ 1 / {xi -yi) \^^ ( {x2-y2) V'^ 

h2\{x- xi) ) V (y - 2/2) / 



(A.13) 



^We take the integral over r' along the real axis. The final answer we obtain by analytic continuation of all 
integrals into complex plane for all variable. 
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In Eq. flA.13|) we found the limit at small values of /12 which contribute to the integral of 
Eq. ()A.5|) . The third integral is equal to 

r ( r'{xi-yi) {''^ ~ '^') {^2 - y2) 



h = / dr' 



00 



(r — r') r' \{x — Xi) {x — r' — yi) ) \{y — I/2) (x — r' — X2) 



hi \ {x-xi) ) I {y-y2) I 



(A. 14) 



Substituting Eq. ()A.11|) and Eq. ()A.13|) into Eq. ()A.6|) we finally obtain the result for the r.h.s. 
of Eq. dSm, namely, 



Q,5 .a+ioo dfi^ I ra+ioo d h2 1 [ {h + h2)\'^ 



J a— too Ani hi \ hih 



(A.15) 



t-2 



' {xi-yi){x-y) \''' I {x2-y2)ix-y) y' ^ ( (xi -yi) ^' / (^2 -^2) ^'' ^ 
{x - xi) {y - yi) ) \{x - X2) {y - y2) ) \{x-xi)) \{y-y2)) 

( {xi -yi){x-y) '^' f {x2 - y2){x - y) V'^ ^ f (xi -^1) ^' / (x2 - ^2) \ 
\{x - xi) {y - yi) ) \{x - X2) {y - y2) ) \{x-xi)) \{y-y2)) 

The integrals over hi and h2 can be evaluated but we postpone this until we work out the action 
of Laplacians (Eq. (j3.3(Jj) ). 

As was noticed in Section 3, the Born amplitude in the form of Eq. ()A.2|) as well as of Eq. ()3.26p 
satisfy the following equation 

A.A,7B^(x,,;x',y') ^ i^^^'iBA{x,y;x',y') (A.16) 

= al (5(2) {x-x') 5(2) {y-y') + 6^'^ {x-y') 6^'^ {y - x')) 

Thus r(2 1) (see Eq. (j3.3Up ) is obtained by applying operator A^ to Eq. ()A.15|) and multi- 
plying by N^/a^. The observation which helps to simplify the calculation is the following: 

where we can neglect terms that are proportional to hf or h"^, as the dominant contribution to 
Eq. ()A.16|) stems from the region of small /I's. Similarly the contribution originating from the 
integrals Ji and Is can be neglected since 

dx dy ^ ^ hi ^ ^ 

Taking into account Eq. ()A.17|) and Eq. ()A.18|) we obtain 

T2^i{ixi,yi) + ix2,y2)^ix,y)) = ^ 177/^^^ ^ ^ ^ 
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^BA{hi;x,y;xi,yi) 7BA (^2; a;, Xa, 1/2) 



{x - vY 

Now we can easily evaluate the remaining integrals over hi and The result can be written in 
the most economic form introducing a new variable: 

R{x.y;x\y') = i'^'^l'i^ (A.20) 
[x-x y [y-y y 

The vertex reads 

Q,3 1/1 

T2^i{{xi,yi) + {x2,y2) ^ {x,y)) = ^^^^^2 ^2 {^x - yY ^ ^' ^i' ^1) + (^-2^) 

2 3 

- hi^ R{x,y]Xi,yi)\nR{x,y]X2,y2) + - \v? R{x,y]Xi,yi)\T? R{x,y]X2,y2) + 

2 1 \ 

+ - In (x, Xi, yi) \v? R (x, X2, ^2) + ^ ^^"^ ^ ^' ^2, 1/2)) 

So far, we evaluated the contribution of the first term (7^]^) of the full Born amplitude of Eq. (jA.ll) . 
Having added the second term we end up with the final expression for the vertex: 

P // X , / \ I w 1 1 R^^^y'i^i^yi) 1 R{x,y;x2,y2) ^ 

^2^i{{xi,yi) + {x2,y2) {x,y)) = In— T^^^t rX 

2{32'KyN^ [x-yy R{x,y;yi, xi) R{x,y;y2, X2) 

(I [ln^i?(x,2/;xi,?/i) + Ini? (x, Xi, Ini? (x, yi, Xi) + In^ i? (x, yi, Xi)] + 

+ I ln(i? (x, y; Xi, yi) R (x, y; yi, Xi)) ln(i? (x, y; X2, y2) R (x, y; y2, X2)) + (A.22) 
+ I [ln^i?(x,?/;X2,?/2) + Ini? (x, X2, 2/2) Ini? (x, 7/2, 3:2) + In'"^ i? (x, y; ?/2, 3^2) 
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